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We study the feasibility of measuring semi-exclusive photon production 7 + p->7 + yat 
HERA. The cross section of photons produced at large transverse momenta, recoiling off an inclusive 
system Y of limited mass, can without photon isolation cuts be simply expressed in terms of hard 
PQCD subprocesses and standard target parton distributions. With the help of event generators we 
identify the kinematic region where quark and gluon fragmentation processes can be neglected. The 
cross section in this semi-exclusive region is large enough to be measured with an upgraded HERA 
luminosity of C = 100 pb _1 . 

The subprocesses of lowest order in a s are suppressed at low recoil masses My , compared to higher 
order gluon exchange, i.e. BFKL contributions. The distinct My-dependence makes it possible to 
determine experimentally the kinematic range where the higher order processes dominate. 
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I. INTRODUCTION 



Semi-exclusive reactions of the type A + B — > C + Y, where particle C emerges with large transverse momentum 
and the mass of the inclusive system Y satisfies Aqcd *C My <C E cm , provide a new tool for probing the structure 
of hadrons jjj. It is analogous to DIS, ep — > eX, in the sense that particles A and C form an effective current which 
probes the target B at large momentum transfer t = (j>a —pc) 2 , producing the inclusive system Y. Semi-exclusive 
cross sections factorize into a calculable short distance interaction times a standard target parton (6) distribution 
fb/B(xs, ~t) according to 



da 



dt dx$ 



(A + B^C + Y) = J2h/ B (xs, 



^da , ., 



Cd) 



(1) 



where 



-t/(My ~ t) 



(2) 



and the momentum transfer —t serves as factorization scale. For large —t only compact configurations of particles 
A and C participate in the reaction, and their re-interactions with the target are suppressed. The final parton(s) d 
merges with the target spectators to form the inclusive system Y (Fig. la). 







FIG. 1. (a) Semi-exclusive photoproduction of photons. The hard process H is calculable in PQCD. (b) Photons produced 
via parton fragmentation are accompanied by comoving hadrons. (c) The resolved photon leaves remnant hadrons in the beam 
direction, (d) The Bethe-Heitler process is not enhanced at low momentum transfer Q 2 between the electrons. 
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We shall discuss the feasibility of measuring semi-exclusive photon photoproduction, -fp — > 7F at HERA. Here 
the incoming and outgoing photons together form an effective current which probes a parton (quark or gluon) in the 
target with resolution —t (Fig. la). This picture is valid provided that the final photon emerges isolated from all 
hadrons in the inclusive system Y, which is ensured by the kinematic requirement W 2 /M Y ^> 1, where W is the total 
center of mass energy, and My is the mass of the hadronic system Y. Scattering on a single parton implies large 
—t, typically of order My. As we shall see, luminosity limitations impose —t <C My at HERA, implying small x s in 
Eq. (|). 

It is not a priori obvious how stringent the kinematic restrictions need to be to ensure semi-exclusive dynamics. At 
moderate W 2 /M Y , photons may be produced via parton fragmentation, following a hard process like 75 — ► gq (Fig. 
lb). Such a process will typically generate hadrons comoving with the photon, implying W 2 /M| ~ 1. Depending on 
the probability that the photon takes most of the momentum of the fragmenting parton such processes can nevertheless 
be significant at moderate values of W 2 /M Y . Based on the PYTHIA g] and LUCIFER [f| Monte Carlo generators 
we conclude that fragmentation processes are insignificant for W 2 /My > 10 and — t > 4 GeV 2 . 

Resolved photon contributions (Fig lc) are analogously suppressed since hadrons in the beam direction must carry 
little energy for W /M Y to be large. This requires the beam parton to carry nearly all the photon energy, and not to 
emit collinear bremsstrahlung prior to its hard scattering. The pointlike photon component transfers the full beam 
energy into the hard process and will thus dominate at large W 2 /M Y and —t. 

We also estimate the contribution of the Bethe- Heftier process (Fig. Id), where the beam photon is virtual and 
the final photon is emitted from the electron. This process is typically detected because of a large angle recoil of 
the electron. We find that the Bethe-Heitler process is in any case negligible for invariant momentum transfers 
Q 2 < 0.1 GeV 2 between the initial and final electrons. 

HERA data [4,|| on photon production at large transverse momentum have previously been compared with QCD 
calculations [0-3 in terms of an isolation cut P-fyJ, which imposes low accompanying transverse energy in a cone 
around the photon direction. Such a cut removes a large fraction of the photons produced by parton fragmentation, but 
allows resolved photon contributions. The isolation procedure requires modelling the non-perturbative fragmentation 
process. This is avoided for the semi-exclusive cross section (Q), which depends only on standard structure functions 
and perturbatively calculable quantities. 



y t y 




g 



FIG. 2. Compton scattering subprocesses: Low order scattering on a quark (a) and on a gluon (b). At high §/ — t, higher 
order amplitudes like (c) dominate, due to gluon exchange in the f-channel. 
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Low order PQCD subprocess amplitudes which contribute to the hard vertex H of Fig. la are shown in Fig. 2. 
It is important to notice that the lowest order contributions of Fig. 2a,b actually are suppressed in the Regge limit 
s = x s W 2 — > oo at fixed (but large) — i, which is relevant for semi-exclusive dynamics. A (dimensionless) amplitude 
involving the exchange of two particles of spins j\ and j% in the i-channel behaves in the Regge limit as ~ (s/ — t) a , 
where a — ji + ji — 1. Hence the contribution of Fig. 2b, with j\ = 32 = 1/2, has a = 0. The amplitude of Fig. 2a 
also has a = 0. On the other hand, the higher order amplitude of Fig. 2c has j\ = ji = 1 an d thus a = 1. At even 
higher orders logarithmic corrections from repeated gluon ladders are expected to build up hard Pomeron exchange 
in the BFKL approximation, as has been studied extensively for the present process jl2]-|l5| . 

It is interesting to explore at what value of s/ — t the higher order diagrams begin to dominate. We find that the 
BFKL amplitude is numerically important compared to the lowest order contributions in the whole semi-exclusive 
range s/ — t > 10. If the BFKL approximation is reliable (or the amplitude of Fig. He dominates) at such moderate 
values of the subprocess energy then the subprocess cross section will be found to grow (or be constant) in this whole 
range. Alternatively, if the lowest order diagrams dominate, the subprocess cross section will decrease in the lower 
range of s/ — t. This would make it smaller than the BFKL estimate at high §/ — t. 



II. SUBPROCESS CROSS SECTIONS 



The lowest order cross section for jq — > jq (Fig. ||a) is |16 



da LO 2-Ke*a 2 ( § -u 



(3) 



with a = e 2 /4ir. For the lowest order 73 — » jg process (Fig. ||b) we have analogously ]T^ , p^ |, 



d<TLO 
dt 



(19 -> 19) = 2 ^ ' 



647TS' 



[\Mi($, t, u)| 2 + \M[(u,t,s)\ 2 + \M[(t 7 s,u)\ 2 + 20] , 



(4) 



where 



t - u . t t 2 +u 2 



|M!(s,i,M)| 2 = 2 + 2— In - + 



s - 1 



\M[(u,t,s)\ 2 = 2 + 2— In — 



, In 2 - 

ir —t 



n \ 2 



In 2 - + tt z 

u 



§ S 2 +t 2 ^ 2 s \ 2 , A _ 2 



4tt z 



s 2 + t 2 , s 



In — 



s-t 



-t u 



(5) 



It is interesting to compare the size of the BFKL cross section to the LO ones given above. The 73 — ► 79 BFKL 
cross section is to leading logarithmic accuracy given by [0-E5I 



cfcr B FKL / \ 81 a 2 a 4 s ( m-^ 



dt 



G(z) = 



16 576t 2 
dv 



7T —t 



with 



tanh(7ri') 



-00 1 + v* 



(" 2 + \Y 



2(ll + 12 I , 2 )exp[0 X (z,)] , 



X (y) = 2^(1) - 2Re 



(6) 



In Fig. ||we show the dimensionless lowest order cross sections s 2 da/dt for ju —> ju (solid line) and jg — > jg (times 
10, dot-dashed line). The small size of the lowest order 73 — > jg cross section means that it will be insignificant also 
in the physical process jN — » jY, even though the gluon distribution is larger than the quark one at low values of 
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x s . The dashed line shows the BFKL cross section (g) for 75 — > jg, based on iterating higher order diagrams like 
that in Fig. |c. We used the fixed value a s = 0.2 indicated by BFKL fits to HERA and Tevatron data The 
BFKL cross section is comparable to the lowest order one already at s/ — t ~ 10, with their ratio growing roughl y a s 
(s/ — i) 2 . It should be stressed, however, that the BFKL analysis is expected to be reliable only for s/ — t > 100 [ |15[ , 
and that there may be sizeable corrections from NLO corrections fl22]. 
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FIG. 3. Lowest order subprocess cross sections s 2 da/dt for ju — » ju (solid line) and 73 — > 73 (times 10, dot-dashed line) for 
10 < s/ — t < 20. The dashed line shows the BFKL cross section extrapolated to this region of moderate subprocess energy, 
assuming a s = 0.2. Three active flavors have been used in both the calculation of the process 75 — > 75 to LO and the BFKL 
contribution. 



III. SEMI-EXCLUSIVE PHOTON ELECTROPRODUCTION AT HERA 



The ep — > e^Y electroproduction cross section is in the semi-exclusive limit Aq C£) <C — t, My <C W 2 related to the 
subprocesses through 

(ep -» e 7 y) = f l/e (y, Q max ) X! fi/p( x " ~ t ^^ i ~* ^ ' ^ 



dydtdx s ^ ' ' 



Here y is the fraction of momentum that the photon, of virtuality < Q 
Williams function f 1 / e (y, <3 max ) i s given in p3j, 



2 

max 



, takes of the electron beam. The Weizsacker- 



f-y/ e (y,Q 



2 

max 



a 



(l + jl-y) 2 ) ^ 



2 

max 



(1-2/) 



_ 2 2 
m^y 



2m 2 e y 



) 2 

J max 



2 2 



(8) 



In our kinematical range the Weizsacker- Williams scale is identical to the upper limit of the momentum transfer 
between the electrons, as the hard scale —t, which characterizes the production process, is much larger than Q 2 max - 
A detailed discussion of the choice of the Weizsacker- Williams scale can be found in Refs. ^^4|. For the parton 
distributions fi/ p (x s , — t) we use the GRV94 LO parameterizations pq ]. 

In this section we first estimate the kinematic region in which the semi-exclusive production mechanism dominates, 
and then evaluate the physical cross section in this region. 
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A. Background from Fragmentation and Decay in Photoproduction 



Large transverse momentum photons can be produced through quark and gluon fragmentation q,g — > 7 + X, 
following standard hard scattering processes such as photon gluon fusion jg — > qq and gluon bremsstrahlung 75 — > qg. 
A second source of contributions are photons from hadronic decays like 7r° — > 77. In the following we call both 
background contributions 'fragmentation'. Such processes typically give hadrons in the photon direction, and are 
thus suppressed in the limit where My <c W 2 . 

We have used the Monte Carlo event generators PYTHIA j| and LUCIFER j| to estimate the range of W 2 /M^ 
in which semi-exclusive production dominates over fragmentation. These event generators include the direct (semi- 
exclusive) 7<7 — > 7(7 cross section only at lowest order. As we observed in the introduction, the higher order contribution 
of Fig. 2c (and its possible BFKL enhancement) actually dominates the lowest order process (Fig. 2a) at high sj — t. 
Hence comparing the fragmentation background to only the LO semi-exclusive cross-section gives a conservative 
estimate of the kinematic region where semi-exclusive dynamics dominates. It should also be kept in mind that the 
fragmentation contribution can be reduced with the help of photon isolation cuts ^,^,^-|TlJ . 

We evaluate the contributions from the direct process 717 — > jq and that from the background 75 — > qq and 
11 — * 19 processes separately. If the final state contains several photons, we choose the one with largest energy E~j M 
in the CM frame of the photon and the nucleon. The invariant mass My of the remaining particles is then given by 
Ml = W 2 - 2WE^ M . 



4 GeV 2 < -t < 10 GeV 2 



10 GeV 2 < -t < 30 GeV 2 
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FIG. 4. Inclusive cross sections da^p — > jY) / 'd(W 2 / 'My) for the photon with the largest E^ M per event at W = 200 GeV. 
The Grange is 4 GeV 2 < -t < 10 GeV 2 (left), and 10 GeV 2 < -t < 30 GeV 2 (right). The triangles show the direct (-yq -> -yq) 
contribution of Eq. (3) and the squares show photons generated by fragmentation. The open squares and triangles are the result 
of the PYTHIA program and the filled squares and triangles are from the LUCIFER run. An additional cut .001 < x s < .7 
has been applied in both figures. Higher order parton showering and resolved photon contributions are not included. 



In Fig. I we show the PYTHIA and LUCIFER results for W = 200 GeV and two ranges of photon transverse 
momentum, 4 GeV 2 < — t < 10 GeV 2 and 10 GeV 2 < — t < 30 GeV 2 . It may be seen that the semi-exclusive process 
11 ~ ¥ 11 starts to dominate for W 2 /A1 2 ~ 10. We have checked that the same conclusion holds for W — 80 GeV. 
Note that for My 3> —t, the subprocess variable s/ — t of Fig. || is in the semi-exclusive limit simply related to 
IU 2 /M y , 
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B. Background from the Bethe-Heitler Process in Electroproduction 



When the incoming photon is virtual the final photon may be radiated off the electron in the Bethe-Heitler (BH) 
process of Fig. | a. We wish to determine the maximum value of Q 2 for which the BH cross section can be neglected 
compared to that of photon emission from the quark, i.e., the Virtual Compton Scattering (VCS) process of Fig. | b. 
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FIG. 5. Sample diagrams for photon production via the Bethe-Heitler (BH) process (a) and Virtual Compton Scattering 
(VCS) (b). 



At finite Q 2 the lowest order ep —> ejY cross section can be written 
da L o(ep — > e-fY) _ 1 



dx s dW 2 dQ 2 dtd<j) 2(4tt) 4 x s W* (s + Q 2 ) 



^2fq/p(Xs,-t)\M q \ 2 , With 



(10) 



\M q \ 2 = \M^ CS \ 2 + \M^ H \ 2 + 2Re (M^ H M^ GS *) 



Here Q 2 = —q 2 = —(I — I') 2 and s = (q + p) 2 = x s W 2 p in the notation of Fig. |j| (f) is the angle between q 'x p 1 and 

I x £' with the three- vectors given in the CM frame of the virtual photon and the proton, while W ep « 300 GeV is 
the ep CM energy. We also define 

si = (l'+ Pl ) 2 , ui = {l- Pl ) 2 (11) 



which can be obtained from s and u by replacing the quark momentum p(p') with the lepton momentum /(/'). 

The squared VCS and BH matrix elements can now be written as 

-4(47ra) 3 4 



\M v q cs \ 2 = 



Q 2 su 
4(47ra) 3 e^ 

UlUl 



(12) 
(13) 



with the common factor F given by 

F = s 2 + [s + t) 2 + uf + (2W 2 q + in) 2 - 2H2W 2 q + ui) - 2t(W 2 q + u t ) + 2Q 2 (2s + t - 3W 2 q - 2ui) + 3Q 4 



(14) 



andW 2 ={l+p) 2 = x s W t 



ep ' 
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2Re (Mf H MX CS *) = 



8(47ra)"e 



3^3 

1 ) r& 



Q 



2 \si ui I ui si 



+t 



-8(1 -pi' -p' + l-p' l'-p) 



la-p+i 1 ■p)-1(i-p' + i' -p') 

u s 



l-p l-p' l'-p l'-p' 



UUi SUi 



USl 



SSi 



+1-P 
+l-p' 
+l'-p 

+l'-p' 



4?- - 3 + A (SI' -p + sil-p' -(* + ui)l' ■ p') 
uui uui 

^ - 3 + ir (-hi -p-*l' -p' + (s + ui)V ■ p) 

SUi SUl 

3 + — (-si ■ p - ud' -p' + (* + §i)l ■ p') 



USl 



USl 



^ - 3 + -^r(ul-p' + uil' -p-(s + si)l ■ p) 

SSi SSi 



(15) 



The results presented in Eqs. ( |12| ) - (15) have previously been computed, e.g., in Ref. |26|], but with our particular 
choice of kinematical variables we obtained more compact expressions. 
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FIG. 6. Dependence of the total ep — > ejY cross section on the maximal momentum transfer Q4ax- The solid line shows the 
contribution of Virtual Compton Scattering (VCS, Fig. 5b) alone, while the dashed line includes also the Bethe-Heitler (BH, 
Fig. 5a) contribution and the interference term. The kinematic region is defined in the text. 



In Fig. || we compare the VCS to the total VCS+BH cross section for ep — > ejY as a function of Q^, the 
maximum momentum transfer between the electrons. The differential cross section (|l0| ) was integrated over the 
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ranges x s G [0.01, 0.7], W € [40 GeV, 160 GeV], </> e [0, 2tt], -t 6 [4 GeV 2 , 10 GeV 2 ], and Q 2 e [m 2 , Q 2 nax ][} 



C. The Semi- Exclusive Cross Section 



We have found that the process ep — > ejY can be used to study semi-exclusive photon production at HERA in the 
kinematic range W 2 /My ~ 10, Q 2 < 0.1 GeV 2 and — t > 4 GeV 2 . In Fig. |?] we show the cross section 

for two ranges of momentum transfer, 4 GeV 2 < —t < 10 GeV 2 and 10 GeV 2 < —t < 30 GeV 2 . The full curve 
shows the contribution of the LO process ofEq. (§) (Fig. |a). The dashed curve shows the contribution of the BFKL 
subprocess of Eq. ( ft) ( Fig. ||c, plus gluon ladder iterations), together with the corresponding jq — ► sea quark 
BFKL contribution 111, 




FIG. 7. The cross section da(ep -> e-yY) / d{W 2 / M$) of Eq. (16) shown for two ranges of t: 4 GeV 2 < -t< 10 GeV 2 (left) 
and 10 GeV 2 <-t< 30 GeV 2 (right). The solid line shows the contribution of the LO process of Eq. (3) (Fig. 2a) . The dashed 
line shows that of the BFKL process of Eq. (17) (Fig. 2c plus gluon ladder iterations). The horizontal dash-dotted line gives 
the 1 event level per unit of W 2 /My, for a HERA luminosity of 100 events/pb. 



*In order to ensure that the electron propagators are far off-shell we also required si, — ui > — tmin = 4 GeV 2 . This cut is. 
however, irrelevant for Q 2 ^ 1 GeV 2 . 







where £ denotes the quark singlet distribution E = Yl q [fq/p + /<?/?>]■ We use a s — 0.2 in the BFKL cross section 
and assume three (four) active flavors in the low (high) |f|-range. The horizontal dash-dotted line indicates the 1 
event level per unit of W 2 /M Y , given a nominal HERA luminosity of 100 events/pb. 

We should emphasize that the BFKL approximation may not be reliable for W 2 /M Y < 100 [[l5), hence the dashed 
curve is an extrapolation. It is nevertheless interesting to observe that this extrapolation dominates the LO qj — > qj 
contribution over the whole semi-exclusive range W 2 /M Y > 10. This conclusion is insensitive to the i-range and also 
to the range of y (not shown) . 



IV. CONCLUSIONS 



Semi-exclusive processes A + B — > C + Y provide a new tool for investigating hadron structure. Effective currents 
formed by the AC system generalize the virtual photon probe familiar from DIS and can carry charge, flavor, baryon 
and other quantum numbers Qj. Before this tool can be put to use, at least two questions need to be answered: 

(i) How stringent limits A.q CD <C —t, M Y <C W 2 must be imposed in order for the semi- exclusive production mechanism 
to dominate? 

(ii) Can the hard AC vertex be reliably computed using PQCD? 

In this paper we studied the process "fp — > jY, which is especially simple in the sense that both particles A 
and C have a point like component. Based on simulations with the PYTHIA and LUCIFER event generators we 
concluded that semi-exclusive dynamics should dominate for W 2 /M Y > 10 and — t > 4 GeV 2 . Photon emission from 
the electron (the Bethe-Heitler process) is insignificant for incoming photon virtualities Q 2 < 0.1 GeV 2 , and can be 
further suppressed with angular cuts. The semi-exclusive cross section should be measurable at HERA, assuming that 
the subprocess cross section a(jq — > jq) is not smaller than its lowest order (LO) PQCD approximation. 

Point (ii) above is non-trivial, since the semi-exclusive kinematics implies a high energy (Regge) limit for the 
subprocess, s/-f > 1. Little is known about the importance of higher order (HO) PQCD corrections in this limit. In 
the process under study the situation is particularly intriguing since the LO subprocess diagrams shown in Fig. 0a, b 
correspond to qq exchange in the ^-channel. At high subenergies s they are therefore suppressed by a factor 1/sr in 
the cross section compared to the O{o i s ) gluon exchange contribution of Fig. ||c. The latter is, on the other hand, 
just the first term in the series of gluon ladder diagrams which is supposed to build up the BFKL Pomeron in this 

iiMi- 



process 



It is not clear from which value of §/ — t HO contributions like Fig. ||c start to dominate the LO processes of 
Fig. ||a,b. The BFKL approximation has been assumed to be relevant for s/ — t > 100 JlSj . Extrapolating the BFKL 
cross section to lower energies we found (Fig. ^) that it would in fact dominate the LO cross section in the whole range 
of semi-exclusive dynamics, s/ — t> 10. The ratio ohoIo~lo behaves approximately like (W 2 /M Y ) 2 , closely reflecting 
the §/ — t dependence of the respective subprocesses, cf. Eq. (^). A measurement of the (W 2 / M Y ) 2 dependence of 
the cross section will thus directly determine the nature of the dominant i-channel exchange. 

We conclude that the large HERA energy in principle allows accessing the semi-exclusive kinematic region, with its 
double hierarchy of large scales. The limiting factor will be the luminosity. If the cross section is approximately given 
by the lowest order contribution in Fig. 7 then only a restricted range of (W 2 / M Y ) 2 can be studied. Since the higher 
order two gluon exchange (BFKL) contributions fall off more slowly with (W 2 / M 2 ) 2 they will eventually dominate. 
If their normalization is even close to that indicated by the BFKL extrapolation of Fig. 7 there should be a rich 
semi-exclusive phenomenology at HERA, not only for photon but also for meson (tt, p, J/ip) production. Further 
work is needed to estimate the feasibility of measuring charge exchange processes and the three gluon (Odderon) 
contribution to ir° production. 
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